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The  goal  of  this  report  was  to  analyse  the  transonic  compressor  tests  available  in  the  open 
literature  and  to  propose  possible  improvements  in  total  loss  and  turning  correlations  From  the 
work  performed  six  conclusions  could  be  drawn  for  the  achievement  of  better  loss  predictions. 
A  comparison  with  results  of  a  two-stage  compressor  showed  that  the  results  of  the  new 
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spanwise  loss  mixing  procedures  must  be  used. 
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PREFACE 


This  report  is  the  outcome  of  AGARD  Support  Project  T15.  The  support  programme  of  AGARD  to  Southern  Flank 
Nations  started  in  198 1,  its  purpose  being  to  set-up  cooperation  among  member  nations  and  thereby  improve  the  capabilities 
of  Southern  Flank  Countries  in  the  field  of  aerospace  research  and  development.  The  programme  is  intended  to  provide 
assistance  to  member  nations  for  the  purpose  of  increasing  their  scientific  and  technical  potential. 

Project  T15  was  initiated  from  the  work  of  Working  Group  12  of  the  AGARD  Propulsion  and  Energetics  Panel  (PEP). 
The  work  was  mainly  conducted  in  the  Mechanical  Engineering  Department  of  the  Middle  East  Technical  University 
(METU),  with  the  technical  support  of  Vrije  Universiteit  Brussel,  Belgium  and  Iowa  State  University,  USA. 

The  authors  would  like  to  thank  the  National  Delegates  of  the  countries  involved,  the  members  of  the  Propulsion  and 
Energetics  Panel,  the  Turkish  Aircraft  Industries  (TUSA§)  for  their  support  in  the  initial  phase  of  the  project,  and  the  Turkish 
Scientific  and  Technical  Research  Council.  The  members  of  the  Fluid  Mechanics  Group  (METU,  Mech.  Eng.  Dept.)  who  took 
part  in  this  project  on  various  occasions  are  acknowledged. 
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Le  present  rapport  s’inscrit  dans  la  suite  logique  du  Projet  de  Support  T15  de  PAGARD.  Le  programme  de  support 
AGARD  aux  Nations  du  Flanc  Sud  date  de  1 98 1 .  Le  but  du  programme  est  de  promouvoir  une  cooperation  parmi  les  nations 
membres  de  TOTAN  et  d’ameliorer  les  potentialites  des  Pays  du  Flanc  Sud  dans  le  domaine  de  la  recherche  et  des  realisations 
aerospatiales.  Le  programme  doit  servir  d’aide  aux  nations  membres  de  POTAN,  et  leur  permettre  d’accroitre  leurs  capacites 
scientifiques  et  techniques. 
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Vrije  Universiteit,  Bruxelles,  Belgique  et  du  Iowa  State  University,  USA. 
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APPLICATION  OF  MODIFIED  LOSS  AND  DEVIATION  CORRELATIONS 
TO  TRANSONIC  AXIAL  COMPRESSORS 

by 

M.^etin,  A.§,Uger,  Ch.Hirsch  and  G.K.Serovy 


1.  INTRODUCTION 

Throughflow  calculation  programs  are  probably  the  most  important  tool  of  the  compressor  aerodynamic  designer.  The 
main  objective  of  a  throughflow  calculation  is  to  provide  a  spanwise  prediction  of  thermodynamic  and  other  flow  variables  so 
that  suitable  blade  profiles  can  be  selected  to  cope  with  the  variations  of  inlet  angle,  turning,  Mach  number,  etc.  Throughflow 
programs  may  be  either  of  analysis  or  design  type.  In  the  former,  which  is  more  extensively  used,  the  blade  geometry  is  specified 
and  solutions  are  sought  for  the  resulting  flow  pattern.  The  design  method,  which  is  less  commonly  available,  requires  enthalpy 
variations  and  obtains  solutions  for  meridional  velocity  and  hence  flow  angles.  To  determine  realistic  solutions  with  both 
techniques  it  is  necessary  to  provide  reasonably  well  predicted  entropy  gradients  resulting  from  viscous  effects.  Similarly,  the 
solution  obtained  from  an  analysis  program  is  very  dependent  upon  the  predicted  values  of  blade  exit  flow  angles.  Once  loss 
and  deviation  models  have  been  incorporated  into  a  program  it  can  be  used  to  predict  overall  machine  performance. 

Continuation  of  improvements  in  throughflow  analysis  programs  and  in  the  empirical  correlations  necessary  for  them  has 
become  one  of  the  interests  of  the  AGARD  Propulsion  and  Energetics  Panel.  The  Panel  organized  a  Specialist  Meeting  in  1 977 
on  “Throughflow  Calculations  in  Axial  Turbomachinery”.  A  Working  Group  (WG  12)  was  subsequently  formed  on 
“Throughflow  Calculations  in  Turbomachinery”.  This  group  started  work  in  January  1978  and  published  its  report  in  1 98 1 .  A 
part  of  this  report  is  devoted  to  the  work  of  the  WG  12  compressor  subgroup.  The  task  of  the  compressor  subgroup  was  to 
review  many  aspects  of  compressor  correlations.  As  a  result  of  this  work  it  was  concluded  that  no  set  of  correlations  has  general 
validity,  and  this  was  especially  emphasized  for  the  validity  of  off-design  correlations.  This  conclusion  was  reached  for  both  loss 
and  deviation  correlation.  As  general  recommendations  of  the  Working  Group  it  was  stated  that: 

(a)  Care  should  be  taken  to  define  a  constant  set  of  relations  within  the  correlations  in  order  to  avoid  the  duplication  of  effects 
and  influence  of  parameters. 

(b)  Notwithstanding  the  continuous  development  of  three  dimensional  viscous  flow  calculations  and  the  increase  of 
computing  power  of  digital  computers,  it  is  not  likely  that  through  flow  calculations  in  multistage  machines  will  be 
replaced  by  full  viscous  calculations  in  the  near  future.  Therefore  the  need  for  reliable  correlations  will  continue  to  be  of 
great  importance. 

(c)  The  general  trends  with  regard  to  the  influence  of  parameters  (such  as  Mach  number)  are  of  particular  interest.  A 
particular  effort  in  this  direction  should  be  made,  especially  for  off-design  cases. 

It  was  with  the  above  conclusions  and  recommendations  that  further  work  on  loss  and  deviation  correlations  was 
proposed.  The  purpose  was  to  use  transonic  compressor  test  results  of  the  1970s  for  reassessing  the  loss  and  deviation 
correlations.  Comparison  of  the  experimental  results  with  the  existing  correlations  was  foreseen.  It  was  expected  either  to 
correct  or  modify  the  existing  correlations  or  to  develop  new  correlations  for  the  off-design  prediction. 

As  the  first  step  of  the  present  investigation,  existing  loss  correlations  were  reviewed  and  put  into  table  form,  which  it  is 
hoped  will  be  useful  for  users.  In  the  table  a  consistent  set  of  notation  is  used  and  the  relevant  tables  and  figures  are  listed  with 
their  number  as  stated  in  the  reference  from  which  the  correlation  is  adopted.  The  table  is  given  in  Appendix  11  with  the 
notation  used  in  the  table  of  Appendix  I,  together  with  cascade  terminology.  Appendix  II  lists  the  equations  used  in  different 
correlations  together  with  the  calculation  procedure.  The  range  of  applicability  of  each  correlation  is  also  stated  in  the  remarks 
section. 

Suitable  data  reported  in  a  number  of  NASA  and  Pratt  &  Whitney  compressor  tests  were  analysed  to  obtain  consistent 
trends  which  might  result  in  new  correlations.  During  this  process  an  off-design  loss  correlation  was  obtained  for  transonic 
upstream  Mach  numbers.  Design  incidence  correlation  by  NASA  for  plane  cascades  was  modified  for  transonic  upstream 
conditions  and  design  deviation  angle  prediction  of  Carter’s  was  corrected.  More  successful  comparisons  were  obtained  from 
modified  and  new  correlations  when  their  results  were  compared  with  those  of  the  other  loss  and  deviation  methods.  No 
consistent  experimental  data  was  found  at  the  transonic  upstream  conditions  for  off-design  deviation  angle  assessment.  Hence 
several  off-design  deviation  correlations  were  examined  and  it  is  concluded  that  the  correlation  given  by  Creveling  is  most 
successful. 

Section  3  of  this  report  discusses  the  methodology  used  in  data  analysis  and  comparison.  The  set  of  correlations  which  was 
set  up  during  the  course  of  this  study  is  programmed  and  fed  into  a  finite-element  throughflow  code.  The  correlations  were 
tested  by  predicting  the  details  of  the  spanwise  variation  of  total  loss,  deviation  angle,  axial  velocity  and  density.  The  test  case 
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used  for  this  purpose  was  a  two-stage  fan.  A  listing  of  the  loss  and  deviation  subroutine  which  may  be  adapted  to  other 
throughflow  programs  is  given  in  Appendix  IV. 

It  must  be  noted  that  in  all  multi-stage  performance  prediction  calculations  the  computation  goes  off-design  if  an 
inaccurate  prediction  of  deviation  for  an  upstream  blade  row  is  made.  Therefore,  a  computer  code  should  calculate  design  and 
off-design  loss  and  deviations  simultaneously  at  all  times. 
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2.  TEST  DATA  USED  FOR  CORRELATION  STUDIES 

2.1  General 

The  search  for  test  data  of  transonic  compressor  stages  which  might  be  used  for  loss  and  deviation  assessment  work  revealed 
eight  useful  data  sets  (references  [  1  ]  to  [8]).  These  tests  were  carried  out  in  the  test  facilities  of  the  NASA  Lewis  Research  Center 
and  the  Pratt  &  Whitney  Aircraft  Company.  Test  reports  include  aerodynamic  design  parameters,  blade  element  geometry,  and 
blade  element  performance  as  well  as  overall  performance  of  blade  rows.  Some  of  the  necessary  blade  element  data  was  lacking  in 
the  test  reports  prepared  by  Pratt  &  Whitney  Aircraft  Company.  The  first  five  test  reports  were  the  results  of  a  research  program 
on  axial  flow  fans  conducted  by  NASA  Lewis  Research  Center.  All  data  which  was  used  in  this  work  was  obtained  from 
compressor  stages  designed  using  1 970’s  technology.  The  annulus  areas  of  these  compressors  varied  over  a  large  range,  with  tip 
diameters  ranging  between  0.25  m  and  0.4 1  mand  hub  diameters  ranging  between  0. 1 25  m  and  0. 1 5  m.  Blade  elements  used  were 
either  Multiple-Circular-Arc  (MCA)  or  Double-Circular-Arc  (DCA)  type.  All  rotors  were  equipped  with  part-span  vibration 
dampers.  For  more  information  the  reader  can  refer  to  references  [1]  to  [8]. 

2.2  Test  Procedure  and  Data  Uncertainty 

Since  the  main  objective  of  the  work  was  to  investigate  off-design  behaviour,  blade  element  performance  data  from  stall  to 
choke  conditions  should  be  available. 


The  details  of  the  test  facilities  are  given  in  references  [1]  to  [8].  In  all  tests  from  which  the  experimental  data  had  been 
obtained  and  used  in  this  investigation,  the  following  test  procedure  was  employed. 

The  compressor  rotor  was  set  to  a  predetermined  rotational  speed,  which  corresponds  to  the  required  blade  velocity  Mach 
number  component.  The  sleeve  valve  located  in  the  collector  at  the  downstream  of  the  compressor  stage  was  then  adjusted.  Each 
flow  condition  created  by  closing  the  sleeve  valve  (increase  of  back  pressure)  established  a  new  flow  field  within  the  compressor, 
and  hence  it  was  possible  to  obtain  incidence  angle  changes  with  the  associated  changes  in  total  pressure  loss  coefficient  and 
deviation  angle.  This  procedure  was  repeated  in  the  unchoked  region  until  the  stall  point  was  reached.  The  sleeve  valve 
adjustment  was  repeated  for  several  rotational  speed  settings.  Passage  survey  measurements  were  made  at  9, 1 1  or  16  radial 
positions,  depending  on  the  test,  ranging  from  5%  ro  95%  span,  for  a  number  of  rotational  speeds  and  back  pressures.  Data  were 
recorded  downstream  of  rotor  and  stator  blade  rows  and  upstream  from  the  stage. 

At  each  radial  position  total  pressure,  static  pressure,  total  temperature,  and  flow  angle  were  measured.  At  the  downstream 
of  the  stator  circumferential  measurements  of  total  pressure,  temperature  and  flow  angle  were  performed.  Mass  and  energy 
averaging  procedures  were  applied  to  these  circumferential  measurements.  Axial  and  radial  velocities  and  flow  angles  were 
obtained  from  the  mass-averaged  properties.  Data  recorded  at  the  measuring  stations  were  modified  to  give  blade  leading  and 
trailing  edge  conditions  using  theoretical  considerations.  More  information  on  the  analysis  procedures  can  be  obtained  from  the 
relevant  references  [1]  to  [8]. 


The  calculation  of  total  loss  coefficient  in  references  [1]  to  [8]  is  accomplished  through  the  equation  below 
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(Po2)id  is  equal  to  Pq^  in  the  case  of  stator  blade  rows,  and  total  loss  coefficient  can  readily  be  calculated  from  the  averaged  values  of 
the  measured  data.  However,  for  rotors  (Po2)id  is  the  summation  of  Pyj  and  the  stagnation  pressure  rise  due  to  the  radius  change 
along  the  streamline.  Thus  equation  2.1  is  modified  to 
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Equation  2.3  implies  that  rothalpy  remains  constant  along  a  streamline.  In  calculating  the  rotor  total  loss,  equation  2.2  is 
used  and  it  is  assumed  that  the  real  streamlines  coincide  with  th£  design  streamlines.  Measured  properties  P^i ,  P02 ,  T(,i  and  co 
together  with  the  approximated  value  of  r2  are  used  to  calculate  coy.  The  error  made  in  the  calculation  of  coy  largely  depends  on 
how  well  the  value  of  r2  is  approximated.  The  effect  of  the  measurement  uncertainties  on  Pqi  ,  P02 ,  Toi » which  are  tabulated  in 
Appendix  II,  are  used  and  propagated  through  equation  2.2  for  calculating  the  amount  of  percentage  uncertainty  on  coy.  In  this 
calculation  it  is  assumed  that  the  streamline  is  exactly  passing  through  the  radii  tabulated  in  the  test  reports.  The  uncertainty  in  the 
radii  is  estimated  as  ±  0.33  mm  from  reference  [40].  As  the  next  step  the  uncertainty  in  the  value  of  r2  is  sought.  It  is  concluded  that 
including  the  bias  error  correspondmg  to  ±  0.5  mm  probe  mispositioning  relative  to  the  actual  streamline  a  total  uncertainty 
interval  of  ±  1.86%  is  expected  on  ooy.  Although  streamline  mispositioning  error  should  be  more  of  a  bias  type  this  cannot  be 
detected  in  the  experimental  minimum  total  losses,  since  the  minimum  total  loss  values  do  not  show  a  distinct  bias  when  compared 
withKochand  Smithresults(Fig.3.10).Asfortheoff-design_total  loss  variations  a  consistent  expected  trend  is  detected  with  angle 
of  incidence.  Thus  it  is  concluded  that  the  error  imposed  on  coy  due  to  streamline  mislocation  is  not  of  serious  order.  The  details  of 
the  calculations  are  given  in  Appendix  HI.  Due  to  the  lack  of  data  used  in  averaging  procedures,  the  uncertainty  analysis  for  data 
reduction  was  not  possible.  The  uncertainties  in  the  angle  measurements  are  tabulated  in  Appendix  III. 


4 


3.  METHODOLOGY  OF  DATA  ANALYSIS 

3.1  Overall  Loss  Correlation 

Transonic  decelerating  flow  in  the  turbomachine  is  one  of  the  most  complex  flows  which  can  be  found  in  fluid  mechanics. 
This  is  because  a  large  number  of  factors  affect  the  flow  regime  and  therefore  influence  the  magnitude  of  losses.  The  complexity  is 
overwhelming  because  of  the  fact  that  the  effects  of  several  loss  producing  factors  are  not  independent,  but  are  cumulative  and 
interrelated.  A  well  known  example  is  the  shock-boundary  layer  interaction.  In  the  presence  of  such  an  interaction  the  order  of 
magnitude  of  profile  loss  is  closely  related  to  the  strength  of  the  shock  wave  as  well  as  its  location  on  the  blade  surface.  The  losses  at 
the  end  walls  are  due  to  the  interactive  phenomena  between  blade  boundary  layers,  end-wall  boundary  layers,  passage  shock- 
waves,  leakage  flow  and  secondary  flow  effects.  It  may  be  quite  unrealistic  to  single  out  a  loss  source,  find  a  loss  coefficient  for  it 
from  experimental  work,  and  then  add  these  linearly  for  sections  along  the  blade  span  where  the  particular  loss  under 
consideration  is  of  importance.  The  equations  which  govern  the  entropy  accumulation  are  highly  nonlinear,  thus  loss  assessment 
using  superposition  is  not  likely  to  be  realistic. 

In  the  present  study,  due  to  the  difficulties  arising  from  accurate  assessment  of  interaction  of  loss  factors,  loss  sources  are  not 
separately  identified.  Instead,  the  data  is  analysed  so  as  to  determine  whether  there  are  some  meaningful  trends  in  the  variation  of 
total  loss  coefficients  with  flow  and  cascade  geometry  parameters.  Losses  are  considered  as  a  whole  and  expressed  by  a  total 
pressure  loss  coefficient  which  accounts  for  all  types  of  losses  except  the  blockage  effects  of  end  wall  boundary  layers  and 
clearance  flows.  This  approach  is  advantageous  for  performance  estimation  because  the  combined  effect  of  profile,  shock  and 
secondary  flow  losses  are  accounted  for  simultaneously  whereas  assessment  of  loss  sources  separately  leads  to  a  superposition 
type  of  synthesis  which  is  in  many  cases  unrealistic  and  imposes  some  unverified  assumptions. 

In  this  work  total  losses  are  considered  in  two  parts.  The  first  part  is  the  minimum  loss  produced  in  the  blade  row,  which 
corresponds  to  the  minimum  loss  incidence;  the  second  part  being  the  loss  increase  when  angle  of  incidence  is  not  equal  to  the 
value  of  minimum  loss  incidence. 

3.2  Variation  of  Total  Loss  with  Angle  of  Incidence 

The  first  step  was  to  plot  the  data  as  the  total  loss  versus  incidence,  and  seek  a  variation  of  total  loss  with  incidence  which  has  a 
detectable  minimum  near  the  zero  incidence.  This  kind  of  curve,  which  is  referred  to  as  loss  characteristic  of  a  blade  section,  can 
only  be  obtained  if  the  incidence  is  changed  over  a  wide  range,  so  that  a  meaningful  variation  of  total  loss  can  be  detected.  A 
variation  in  angle  of  incidence  can  only  be  obtained  if  the  test  compressor  back  pressure  is  changed  at  constant  speed  .The  number 
of  tests  that  are  available  in  the  open  literature  having  data  suitable  for  off-design  analysis  is  limited.  Although  a  variation  of  weight 
flow  is  expected  with  a  change  of  back  pressure  at  unchoked  flow,  this  variation  is  quite  small.  Thus  it  is  possible  to  obtain  different 
flow  regimes  by  changing  the  rotational  speed  and  back  pressure  independently.  The  plots  of  total  loss  with  angle  of  incidence  for 
stator  blade  rows  exhibit  results  similar  to  the  sample  given  in  figure  3.1.  As  it  is  clearly  seen  it  is  rather  difficult  to  indicate  the 
minimum  loss  point  from  such  data.  The  test  reports  also  mention  as  a  conclusion  that  for  stator  the  minimum  loss  range  extends 
8°  around  the  minimum  loss  point  which  is  rather  a  large  range.  Thus  the  flatness  of  the  loss  characteristic  around  the  minimum 
loss  directly  affects  the  process  of  determining  the  minimum  loss  incidence.  Therefore  all  stator  data  were  excluded  from  the 
analysis  and  rotor  data  were  used  for_die  basis  of  design  and  off-design  total  loss  analysis.  Some  of  the  loss  characteristics  of  the 
rotors  also  exhibited  variations  of  coj  with  i  which  does  not  allow  definition  of  a  rational  minimum  loss  point.  In  these 
circumstances  the  data  were  not  used. 

Six  sample  loss  characteristic  are  given  in  figures  3.2  to  3.7.  The  samples  are  chosen  from  three  different  spanwise  locations. 
The  first  three  samples  are  for  DCA  blading  and  the  remainder  are  for  MCA  blading.  The  blade  element  section  parameters  and 
cascade  geometry  are  given  in  each  figure.  There  was  difficulty  in  finding  the  variation  of  total  loss  with  incidence  for  a  given  blade 
and  cascade  geometry  at  a  single  constant  inlet  mach  number.  Therefore,  as  is  seen  in  the  sample  figures,  inlet  Mach  numbers 
corresponding  to  the  data  points  are  presented  as  intervals.  The  Mach  number  range  for  each  curve  is  narrow  enough  to  assign  an 
arithmetic  mean  Mach  number  value  for  each  curve. 

It  was  a  formidable  task  requiring  careful  observation  of  the  curve  trend  and  engineering  intuition  to  fit  curves  to  the  data 
points  and  fix  the  minimum  loss  point  of  each  curve.  This  process  gave  the  minimum  loss  incidence  i*  p  and  co*  p  values.  Since  the 
main  objective  was  to  find  off-design  correlation  for  transonic  blade  sections,  it  was  decided  to  check  whether  both  experimental 
minimum  loss  and  experimental  minimum  loss  incidence  can  or  can  not  be  determined  using  existing  correlations. 

3.3  Minimum  Loss  Incidence  Prediction 

It  is  well  known  and  can  be  seen  very  clearly  from  figures  3 .2  to  3 .7  that  minimum  loss  incidence  angle  prediction  is  one  of  the 
most  important  issues  for  the  assessment  of  off-design  performance.  It  is  essential  to  predict  an  accurate  value  for  this  incidence. 

The  incidence  angle  can  be  measured  from  blade  suction  surface  or  from  mean  camber  line  (see  Appendix  I).  It  was  found 
convenient  to  use  incidence  angle  definition  referenced  to  the  mean  camber  line,  in  accordance  with  NASA  incidence  angle 
correlations.  However,  it  is  always  possible  to  convert  the  angle  of  incidence  based  on  suction  surface  i^^  into  the  one  based  on 
mean  camber  line  i,  if  the  geometry  of  the  leading  edge  is  known. 

The  minimum  loss  design  incidence  angle  correlation  presented  in  NASA  SP-3  6  [9]  covers  a  wide  range  of  parameters  which 
were  carefully  examined  for  subsonic  flow  regimes  in  two  dimensional  cascades.  NASA  minimum  loss  incidence  angle 
correlations  also  cover  corrections  for  three-dimensional  rotor  cascades  against  relative  inlet  Mach  number  for  DCA  and  NAC  A 
65-(A10)  series  blade  profiles. 
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Since  the  compressor  test  data  used  in  this  investigation  included  3-D  effects  it  was  decided  to  compare  the  NASA- 2D 
prediction  with  the  experimentally  determined  i*p  values.  Figure  3.8  shows  the  variation  of  i*  p  -  I^d-nasa  with  inlet  Mach 
number.  The  figure  shows  that  i*  p  -  i^D-NASAis  less  for  MCA  blades  than  for  DCA  blades.  The  other  observation  from  the  figure 
is  that,  for  both  types  of  blade  profiles,  the  difference  increases  with  inlet  Mach  number.  Straight  lines  are  fitted  separately  through 
the  data  points  of  DCA  and  MCA  blades.  The  following  equations  are  obtained  for  correcting  the  NASA-2D  minimum  loss 
incidence  prediction. 

For  DCA  blade  profiles 

iSoR  -  i2D  =  0.7238M,  +  7.5481  (3.1) 

For  MCA  blade  profiles 

icoR  -  1*0  =  1.3026M,  +  5.7380  (3.2) 

where  iJoR  in  the  above  equations  is  the  corrected  value  of  minimum  loss  incidence.  Figure  3.9  shows  a  comparison  of  NASA-3D 
minimum  loss  incidence  angle  prediction  method  and  the  presently  modified  NASA-2D  method.  A  random  sampling  of  a 
number  of  DCA  profiles  operating  at  different  flow  conditions  have  shown  that  the  NASA  3-D  prediction  is  close  to  the  present 
prediction  around  inlet  Mach  number  1.0.  However,  at  low  and  high  Mach  numbers,  the  discrepancy  between  the  present 
prediction  method  and  NASA-3D  predictions  increases  to  as  much  as  60%.  The  correction  applied  to  NASA-2D  correlation  in 
reference  [9]  for  3-D  effects  and  inlet  Mach  number  exhibits  no  Mach  number  influence  between  inlet  Mach  number  0.25  and 
0.55.  The  correction  then  increased  with  inlet  Mach  number.  This  increase  continues  approximately  between  inlet  Mach 
numbers  0.55  and  0.95.  The  3-D  cascades  used  to  derive  the  NAS  A-3D  correlation  were  probably  choked  at  inlet  Mach  number 
0.95.  Thus,  a  constant  value  of  correction  is  indicated.  It  can  therefore  be  argued  that  the  variation  of  correction  between  inlet 
Mach  numbers  ~  0.6  and  ~  1.0  is  not  the  same  as  it  is  established  in  this  work.  This  leads  to  a  discrepancy  ini*  predictions  for  M^ 
values  less  than  1 .0  in  Figure  3.9.  For  inlet  Mach  numbers  greater  than  1 .0,  the  NASA  correction  is  no  longer  valid.  In  the  present 
analysis  3-D  and  Mach  number  effects  cannot  be  separated. 

3.4  Design  Total  Loss  Prediction 

As  it  is  mentioned  above,  the  minimum  loss  points  of  the  experimental  loss  characteristics  are  fixed  by  fitting  curves  to  the 
available  test  points  present  at  an  approximately  constant  inlet  Mach  number.  These  fixed  minimum  loss  values  are  then 
compared  with  the  existing  design  loss  prediction  methods  which  include  corrections  for  transonic  effects.  The  comparison 
between  experimental  and  calculated  total  losses  using  different  prediction  methods  can  be  seen  in  Figure  3.10.  Four  different 
design  loss  prediction  methods  are  used.  Koch  and  Smith  [17],  Swan  [16],Dettmering  [34],  Jansen  and  Moffatt  [18].  Anumberof 
calculations  have  suggested  that  the  method  proposed  by  Koch  and  Smith  is  the  most  satisfactory  one  in  determining  the  design 
loss  for  transonic  cascades.  This  method  in  fact  is  the  most  complete  one,  accounting  for  many  parameters  such  as  stream  tube 
contraction,  blade  surface  roughness,  etc.  The  method  uses  an  equivalent  sand  roughness  term  in  the  calculation  of  profile  losses 
which  affects  the  profile  loss  magnitude.  It  is  concluded  in  [1 7]  that  an  accurate  assessment  of  profile  loss  can  only  be  obtained  by 
the  evaluation  of  surface  roughness  with  utmost  care. 

3.5  Off-Design  Total  Loss  Prediction 

The  work  described  above  on  the  prediction  of  minimum  total  loss  w*  and  minimum  loss  incident  i*  showed  that  the 
experimentally  obtained  values  of  these  quantities  can  be  determined  with  less  error  using  the  correct  form  of  NASA-2D 
correlation  for  minimum  loss  incidence  angle  Koch  and  Smith  correlation  for  minimum  total  loss. 

The  experimental  loss  characteristics  were  then  replotted  as  (wt  -  (o^ )  v.  (i  -  i*),  keeping  in  mind  that  i*  and  (o^  can  be 
determined  from  the  a^ve  r^ntioned  correlations.  Initially  the  loss  characteristics  for  MCA  and  DCA  blade  profiles  were 
replotted  on  the  same  (wj  —  coj)  v.  (i  —  i*)  plane.  It  was  then  realised  that  there  are  observable  differences  in  the  slopes  of  the 
characteristics  for  MCA  and  DCAblades.  It  was  therefore  concluded  that  the  magnitude  of  losses  depends  on  the  profile  shape. 
The  plot^are  r^eated  separately  in  Figures  3.11  and  3.12  for  MCAand  DCA  blade  profiles  respectively.  The  origin  in  these 
figures,(a)T  “  (^>7)  =  0  and  (i  -  i  )  =  0,  corresponds  to  experimentally  determined  values  of  minimum  loss  and  minimum  loss 
incidence. 

The  figures  show  that  at  off-design  incidences  the  total  loss  parameter  w*  depends  largely  on  the  inlet  Mach  number. 
Although  casc^e  geometrical  variables  such  as  solidity,  camber  angle,  maximum  thickness,  etc.,  influence  to  a  large  extent  the 
minimum  loss  coj ,  no  systematic  influence  of  the  cascade  geometrical  parameters  on  the  off-design  loss  can  be  detected. 

Acomparison  between  the  off-design  losses  for  DCAand  MCA  blades  showed  that  the  losses  produced  by  DCA  type  of 
blading  are  higher  compared  to  MCA  type  of  blading.  This  is  because,  for  MCA  type  of  blading,  camber  angle  and  blade  thickness 
are  smaller  at  the  enterence  region,  thus  providing  lower  suction  surface  Mach  numbers  just  ahead  of  the  passage  shock  for  a  given 
range  of  Mach  number. 

Acareful  examination  of  the  data  points  shows  a  very  consistent  inlet  Mach  number  dependence  at  the  negative  values  of 
(i  -  i*).  As  the  inlet  Mach  number  increases  at  a  constant  negative  value  of  (i  -  i*)  the  total  off-design  loss  increases.  On  the  right 
hand  side  of  (i  -  i*)  =  0,  Mach  number  dependence  is  not  very  distinct  because  of  the  large  scatter  of  data  points.  However,  the 
increase  of  off-design  loss  with  inlet  Mach  number  can  still  be  detected  if  least  square  fits  are  applied  to  data  points. 
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It  should  once  more  be  strongly  emphasised  that  the  experimental  data  points  in  Figures  3.11  and  3. 1 2  are  for  MCA  and 
DCA  airfoil  sections  with  a  large  range  of  cascade  geometries.  The  ranges  of  solidity,  camber  angle,  stagger  angle,  and  maximum 
thickness  to  chord  ratio  are  indicated  in  the  figures.  This  indeed  shows  that  the  inlet  Mach  number  and  the  type  of  blading  are  the 
two  strong  governing  factors  in  determining  off-design  losses  at  transonic  speeds. 

Datafor{i  -  i*)  <  Oand(i  -  i^)  >  0_are  treated  separately.  The  available  data  for  each  inlet  Mach  number  are  used  for  least- 
squares  fit  to  a  function  in  the  form  (Ot  -  (0*  =  Cn,(i  -  i*)".  A  large  number  of  trials  revealed  thm  curves  with  exponent  n  equal  to 
2  give  the  best  fit  to  the  data  points  for  all  inlet  Mach  numbers  at  both  sides  of  the  origin.  Thus  coj  -  co*  =  Cn,(i  -  i*f  is  chosen 
and  values  of  c^,  for  all  inlet  Mach  numbers  are  determined.  The  fitted  curves  are  shown  in  Figures  3.13  and  3.14  for  multiple 
circular  arc  and  double  circular  arc  blades  respectively. 

The  relations  between  inlet  Mach  number  and  the  coefficient  c^  are  given  by  Figures  3.15  and  3.16.  Figure  3.15  shows  the 
variation  of  c^^  with  inlet  Mach  number  for  MCA  blade  profiles.  The  variation  can  be  approximated  by  a  straight  line.  A  straight 
line  is  fitted  to  the  data  points  using  least  square  curve  fitting  technique.  The  slope  ofthe  straight  line  for  (i  -  i*)  <  0  is  greater  than 
that  of  (i  -  i*)  >  0.  This  means  that  for  negative  incidences  the  total  loss  increases  faster  with  Mach  number  at  off-design.  The 
same  trend  can  also  be  seen  from  F igure  3. 1 6  for  DCA  blade  profiles.  Double  circular  arc  profiles  exhibit  higher  off-design  losses 
at  the  same  inlet  Mach  numbers  compared  to  the  MCA  blades.  Variation  of  c^  with  inlet  Mach  number  is  given  by  the  equations 
given  in  Table  1 . 

The  equations  in  Table  together  with  the  modified  NASA-2D  correlation  for  minimum  loss  incidence  and  Koch  &  Smith 
design  loss  correlation  may  be  used  to  predict  the  off-design  loss  performance  of  blade  sections  for  design  and  performance 
prediction  purposes.  For  calculating  the  total  pressure  loss  coefficient  for  a  given  blade  and  cascade  geometry,  inlet  flow  direction 
and  inlet  Mach  number  one  can  proceed  as  follows; 

Table  1  Variation  of  Coefficient  Cj^  with  inlet  Mach  Number 
Base  equation  (Oj  -  (o*  =  c^{i  -  i*)^ ,  Mi  >  0.5 


Type  of 
Blade  Profile 

(i  -  i*) 

Equation 

MCA 

<  0 

c„  =  0.02845M,  -  0.01741 

>  0 

c^  =  0.00363Mi  -  0,00065 

DCA 

<  0 

=  0,05336Mi  -  0,02937 

>  0 

c„  =  0.00500M,  -  0.00075 

1 .  Calculate  the  minimum  (design)  loss  coefficient  (o*  using  Koch  &  Smith  [17]  correlation. 

2.  Calculate  minimum  loss  incidence  i*  using  modified  NASA-2 D  correlation  (see  Section  3.3 ) 

3.  Calculate  the  value  of  coefficient  c^,  from  the  equations  given  in  Table  1 . 

4.  Calculate  off-design  loss-coefficient  from 

0)t  =  CO*  +  cji  -  i*)^ . 


Equation  3.3  can  also  be  used  to  estimate  stall  and  choke  angle  of  incidences  if  they  are  taken  as  the  incidences,  where  the 
losses  are  twice  their  minimum  values.  Equation  3.3  can  now  be  written  in  the  form 


(3.4) 


where  c^^  can  be  obtained  from  Table  1 .  The  value  of  c^^^  for  predicting  stall  incidence  is  obtained  from  equations  given  for 
(i  —  i*)  >  0,  whereas  for  choke  incidence,  equations  for  (i  —  i*)  <  0  are  used.  The  values  of  i*  and  co*  are  calculated  from 
modified  NASA-2D  and  Koch  &  Smith  correlations. 


3.6  Verification  of  the  Total  Loss  Correlation  Set 

The  four  off-design  correlations  of  Creveling  et  al.  [  1 9],  Swan  [16],  Howell  [20]  and  Jansen  &  Moffatt  [  1 8]  are  used  to  predict 
the  off-design  loss  and  the  results  are  compared  with  the  measured  values  of  (coj  -  co*).  For  these  calculations  samples  are  taken 
arbitrarily  from  different  blade  profiles,  cascade  geometry  and  upstream  conditions.  Figures  3.17  to  3.20  show  these 
comparisons.  From  the  figures  it  may  be  concluded  that,  although  a  slight  under-prediction  exists,  the  best  prediction  is  made  by 
Swan’s  off-design  correlation.  The  large  over  predicted  values  by  Howell’s  correlation  are  attributed  to  the  wrong  assessment  of 
fluid  angles  at  design  and  off-design  operations,  and  to  the  lack  of  any  corrections  for  Mach  number  effects.  The  overpredictions 
seen  in  Jansen  &  Moffatt  off-design  loss  calculations  are  most  probably  due  to  the  errors  in  assessing  the  stall  and  choke  angle  of 
incidences.  Creveling’s  off-design  loss  prediction  method  is  more  successful  compared  to  the  above  two.  The  discrepancies  from 
the  experimental  values  are  attributed  to  the  non-co  ntinuous  variation  off-design  loss  correction  with  inlet  Mach  number.  Swan’s 
correlation  is  more  successful  since  it  accounts  for  the  momentum  thickness  change  with  the  inlet  Mach  number. 
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F igure  3.2 1  shows  a  comparison  of  the  experimental  off-design  loss  and  the  off-design  loss  predicted  by  the  new  correlation. 
It  is  clearly  seen  that  the  new  off-design  correlation  predicts  the  off  design  losses  most  satisfactorily.  The  results  are  condensed 
around  the  one  to  one  correspondence  line  with  a  low  degree  of  scatter. 

Two  other  verification  calculations  are  also  performed  by  using  the  new  loss  correlation  set.  The  loss  coefficient,  including 
design  and  off-design  parts,  is  calculated  for  a  number  of  blade  and  cascade  geometries  and  inlet  conditions.  The  calculated  loss 
coefficients  are  then  compared  with  the  corresponding  experimental  values.  The  first  comparison  is  given  in  Figure  3.22.  This 
figure  is  obtained  from  sample  operating  conditions  selected  arbitrarily  among  the  data  used  for  obtaining  the  new  correlation  set. 
The  figure  shows  the  variation  of  percent  deviation  of  the  calculated  value  from  that  of  the  experimental,  with  inlet  Mach  number. 
It  is  clearly  seen  that  there  is  no  bias  with  regard  to  either  inlet  Mach  number  or  the  shape  of  the  blade.  All  points  fall  inside  ±  30% 
deviation  from  the  experimental  value.  The  range  of  incidence  angles  of  the  data  points  presented  in  Figure  3.22  is  between  -  3 
and  +15  degrees.  The  second  check  for  the  new  correlation  set  was  made  by  utilising  the  unused  data.  Percent  deviation  of 
calculated  loss  coefficient  from  experimental  value  is  plotted  against  inlet  Mach  number  in  Figure  3.23.  A  larger  deviation  is 
expected,  since  this  data  was  not  used  in  the  correlation  assessment  work  due  to  the  large  scatter.  The  data  points  with  higher 
percent  deviations  mostly  belong  to  the  near  stall  operating  conditions.  The  figure  shows  that  about  60%  of  the  data  points  fall 
inside  the  ±  30%  deviation  band. 

3.7  Design  Deviation  Angle  Correlation 

The  data  available  in  references  [1]  to  [8]  were  used  for  deviation  angle  correlation  assessment.  Measured  deviations  were 
plotted  against  angle  of  incidence  for  the  prescribed  inlet  Mach  number  ranges.  A  large  number  of  such  plots  were  made  for 
stators  and  rotors  at  various  spanwise  locations  [40].  The  result  of  careful  inspection  of  all  the  plots  showed  that  the  effect  of  angle 
of  incidence  on  deviation  angle  does  not  show  any  consistent  trend.  A  sample  plot  is  shown  in  Figure  3.24.  No  influence  of  the 
Mach  number  was  also  seen.  For  rotor  blade  elements  a  slightly  increasing  trend  of  deviation  angle  with  positive  incidence  is 
detected  in  some  of  the  tests.  However,  it  was  not  possible  to  detect  this  rise  for  stator  blade  sections.  Therefore,  a  single  deviation 
angle  was  calculated  from  the  arithmetic  mean  of  the  data  points  of  different  inlet  Mach  numbers  and  incidence  angles  of  the 
known  cascade  geometry.  It  was  then  decided  to  compare  these  mean  deviation  angles  with  the  deviation  angles  calculated  from 
correlations  [42].  Figure  3.25  shows  the  comparison  of  these  mean  measured  deviation  angles  with  the  deviation  angles  calculated 
from  Carter’s  rule.  The  figure  implies  that  Carter’s  rule  underestimates  the  deviation  angle  to  some  extent,  especially  at  the  higher 
deviation  values.  In  Figure  3.26  a  similar  comparison  for  NASA/N  AC  A  deviation  is  given.  This  correlation  is  better  in  predicting 
the  experimental  deviations  obtained  from  the  compressor  tests  with  a  larger  scatter.  In  Figures  3.25  and  3.26  data  points  for  both 
MCA  and  DCA  blades  exist.  Since  a  more  consistent  trend  is  observed  in  the  discrepancy  between  measured  and  calculated 
values  with  Carter’s  rule,  a  correction  to  the  Carter’s  rule  is  obtained  by  fitting  a  least  square  second  order  curve  to  the  data  points, 
as  follows: 


6*  -  - 1.099379  +  3.01866car  “  0.19886*^  (3.5) 

The  discrepancy  between  6car  and  6* ,  and  therefore  the  required  correction  can  be  attributed  to  the  transonic  and  3-D  effects 
which  can  not  be  separated  in  the  present  analysis.  The  value  of  coefficient  m  in  Carter’s  rule  was  obtained  by  assuming  circular 
arc  mean  camber  line.  The  comparison  of  modified  Carter’s  rule  and  experimental  mean  deviations  is  shown  in  Figure  3.27. 

3.8  Off-Design  Deviation  Angle  Correlation 

Since  no  consistent  off-design  deviation  angle  data  could  be  extracted  from  the  available  compressor  tests,  it  was  decided  to 
use  throughflow  calculations  for  the  assessment  of  the  best  off-design  deviation  correlation.  For  this  purpose,  off-design  deviation 
angle  correlations  proposed  by  Jansen  &  Moffatt  [18],  Swan  [16],  Creveling  [  1 9]  and  Howell  [20]  were  evaluated.  Systematic  test 
runs  were  performed  for  determining  the  most  successful  off-design  deviation  calculation  technique.  As  a  result  of  these 
calculations  Creveling’s  off-design  evaluation  method  is  found  more  successful  than  the  others.  It  must  be  noted  that  some 
convergence  problems  were  experienced  during  the  computations  performed  by  correlations  other  than  Creveling’s  more 
detailed  information  with  regard  to  the  computations  given  in  Section  4. 
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Figure  3.1  Loss  Characteristic  of  a  Sample  Blade  Element  Section 
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Figure  3.2  Variation  of  Total  Losses  with  Incidence 
(42.5  %  Span) 
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Figure  3.3  Variation  of  Total  Losses  with  Incidence 
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Figure  3.4  Variation  of  Total  Losses  with  Incidence 
(50  %  Span) 
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Figure  3.5  The  Variation  of  Total  Losses  with  Incidence 
(10  %  Span) 
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Figure  3.6  The  Variation  of  Total  Losses  with  Incidence 
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Figure  3.7  The  Variation  of  Total  Losses  with  Incidence 
(40  %  Span) 
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Figure  3.8  Deduced  Variation  of  Transonic  Rotor  Design  Incidence 
Angle  Minus  Low-Speed  Two-Dimensional -Cascade-Rule 
Design  Incidence  Angle  [9]  with  Relative  Inlet  Mach 
Number  For  DCA  and  MCA  Profiles 


100 

80 

60 

40 

20 

0 

20 

40 

60 

80 

100 

0. 


1 

i 

A 

i 

Ik 

1 

A 

^  A 

o 

i 

'  rC 

L 

C 

t  vyO 

qOCLQ 

c 

\  - 

- s 

>  c 

)  ^ 

b  ^ 

O  ^2D-C0R^^exp 
^  ^3D-NASA^^exp 

5  0.6  0.7  0.8  0.9  1.0  1.1  1.2  1.3 


3.9  Comparison  of  Design  Angle  of  Incidence  Calculated  by  NASA-3D  [9]  and  by 
NASA-2D  with  Modification  for  3-D  and  Transonic  Cascades 


17 


EXPERIMENTAL 

Figure  3.10  Comparison  of  Design-loss  Coefficient  Calculated 
by  Various  Correlations  with  Experimental  Results 
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Figure  3.11  Variatioa  of  Off-Desiga  Losses  with  Design  Incidences 
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Figure  3.12  Variation  of  Off-Design  Losses  with  Off-Design  Incidences 
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Figure  3.13  Curve  fit  through  the  Data  Points  of  MCA  Blade  Profiles 


Figure  3.14  Curve  fit  Through  the  Data  Points  of  DCA  Blade  Profiles 
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Figure  3.17  Comparison  of  Creveling's  [19]  Off-Design  Loss 
Correlation  with  Experimental  Results 


T^calculated 


24 


0.45 


0.40 


0.35 


0.30 


0.25 

13 

'  0.20 
h- 

s. 

0.15 


0.10 


0.05 


0.05  0.10  0.15  0.20 

(wj  -  “j)experiiTiental 

Figure  3.18  Comparison  of  Swan's  [16]  Off-Design  Loss 
Correlation  with  Experimental  Results 
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Figure  3.20  Comparison  of  Jansen  &  Moffatt's[18]  Off-Design  Loss 
Correlation  with  Experimental  Results 
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Figure  3.24  A  Sample  of  Deviation  Versus  Angle  of  Incidence 
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Figure  3.25  Comparison  of  Measured  Deviation  Angles 

with  the  Predictions  of  Carter's  Correlation 
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4.  APPLICATION  OF  MODIFIED  LOSS  AND  DEVIATION  CORRELATION  SET  TO  PERFORMANCE 
PREDICTION 


As  a  result  of  the  work  explained  in  the  last  section  a  set  of  correlations  were  selected  which  may  be  used  for  performance 
prediction  of  transonic  compressors.  These  are: 


Minimum  loss  incidence  angle:  Modified  NASA-2D 
Design  total  loss:  Koch  and  Smith 

Off-design  total  loss:  New  Correlation  (Table  1 ) 

Design  deviation  angle:  Modified  Carter’s  rule 

Off-design  deviation  angle:  Creveling 


It  should  be  noted  that  the  design  loss  and  deviation  are  the  values  which  occur  at  the  minimum-loss  incidence  angle. 
Therefore  blade  sections  will  be  operating  at  off-design  condition  if  i*  i.  This  condition  almost  always  occurs  in  performance 
prediction  calculations  because  of  the  iterative  procedures  and  inherent  errors  of  the  computational  method.  Especially  during 
the  convergence  phase  of  the  computation,  the  blade  sections  would  be  operating  far  from  design.  Therefore,  some  stops  must  be 
included  in  computer  codes  in  order  to  control  the  loss  and  deviations  when  the  blade  rows  either  stall  or  choke  prematurely. 


The  throughflow  computer  code  used  for  performance  prediction  calculations  was  developed  in  the  Mechanical 
Engineering  Department  of  the  Middle  East  Technical  University  [37].  The  code  uses  a  finite-element  method  for  determining 
the  flow  on  axisymmetric  hub-to-shroud  surfaces.  The  code  requires  empirical  correlations  for  determining  the  losses  and 
turning  through  the  blade  rows.  Interblade  calculation  stations  are  allowed  and  blade  blockage  and  lean  angle  are  distributed  to 
these  locations  from  the  geometrical  data.  Turning  and  loss  are  also  estimated  for  the  interblade  nodes  in  the  stream  wise  direction. 
In  order  not  to  complicate  the  problem  by  introducing  additional  unknowns,  no  end  wall  boundary  layer  calculation  procedure  is 
used.  Instead,  experimental  values  of  the  end  wall  boundary  layer  blockage  are  introduced.  This  is  done  by  displacing  the  hub  and 
tip  walls. 

A  systematic  investigation  was  performed  on  the  influence  of  the  stream  wise  inter-blade  distribution  of  loss  and  deviation 
on  flow  properties.  Both  linear  and  exponential  distributions  of  loss  and  deviation  were  evaluated.  Different  ways  of  distributing 
the  loss  and  deviation  to  the  inter-blade  stations  had  a  negligible  effect  on  the  calculated  flow  property  distributions.  Exponential 
distribution  with  the  exponent  less  than  one  gave  convergence  problems.  The  exponent  equal  to  two  was  successful  with  results 
almost  the  same  as  the  linear  distribution  i.e.,  exponent  equal  to  one.  A  linear  distribution  of  loss  and  turning  is  preferred  and  used 
in  the  performance  prediction  calculations. 

As  the  test  geometry,  a  NASA  two  stage  fan  is  used.  This  test  case  is  well  documented  in  references  [26]  and  [4 1  ].  The  fan  has 
low  aspect  ratio  blading  particularly  in  the  first  rotor  (aspect  ratio  being  equal  to  1.56).  This  test  case  is  quite  challenging  with 
regard  to  performance  prediction,  especially  off-design  as  illustrated  in  AGARD  WG 1 2  report  [11].  The  test  case  is  one  of  the  few 
well-documented  two-stage  compressor  test  cases,  allowing  comparisons  for  multiple  blade  rows. 


As  the  first  check  of  the  computer  code,  calculations  were  performed  at  70, 80, 90  and  100  percent  design  speeds  using  the 
measured  total  loss  and  deviations  given  in  reference  [26].  The  axial  velocity  and  density  variations  in  the  spanwise  direction  were 
found  to  be  satisfactory  at  all  speeds.  Later  computations  were  performed  using  measured  turning  angles  and  correlations  for  loss 
prediction.  Computations  were  also  performed  using  correlations  for  turning  angle,  and  measured  values  of  losses  as  input.  These 
calculations  showed  that  correct  prediction  of  turning  angle  at  design  and  off-design  is  more  important  than  the  correct  prediction 
of  losses. 


For  selection  of  the  most  suitable  off-design  deviation  correlation,  test  computations  were  performed  using  measured  total 
losses  as  input  with  the  various  off-design  deviations  mentioned  in  Section  3.8.  As  a  result  of  these  computations  Creveling’s  off- 
design  deviation  correlation  was  found  to  be  the  most  satisfactory  one. 

The  experimental  and  predicted  values  given  in  F  igures  4. 1  to  4.4  are  the  results  of  the  computations  performed  using  the  set 
of  loss  correlations  recommended  as  a  result  of  the  work  described  in  this  report.  The  design  and  off-design  correlations  are  used 
simultaneously  in  the  calculations.  The  correlation  sub-program  is  given  in  Appendix  IV.  The  test  point  calculated  is  at  100% 
design  speed.  The  comparisons  of  experimental  and  theoretical  results  of  the  first  stage  rotor,  first  stage  stator,  second  stage  rotor, 
and  second  stage  stator  are  given  in  Figures  4.1, 4.2,  4.3  and  4.4  respectively.  In  each  figure  the  variation  of  axial  velocity  and 
density  with  percent  span  is  plotted  both  u  pstream  and  downstream  of  the  blade  row.  On  the  same  figure  the  variation  of  predicted 
loss  and  deviation  with  percent  span  is  also  displayed.  Percent  span  indicated  in  the  figures  is  measured  from  the  hub.  As  expected, 
examination  of  the  figures  shows  that  the  correlation  set  can  not  be  regarded  as  very  successful  in  predicting  the  details  of  flow  near 
the  hub  and  casing.  Because  no  explicit  spanwise  treatment  of  the  secondary  flow  losses  exists  this  mis-p  rediction  extends  towards 
the  mid-span  especially  in  rotor  blades.  The  discrepancies  are  greater  at  the  tip  of  the  rotor  blades  due  to  the  complicated  flow 
structure  of  leakage,  end  wall  boundary  layer  and  transonic  flow  effects  which  can  not  be  completely  accounted  for  in  the 
predicted  total  loss  coefficients.  However,  except  at  the  end  wall  regions,  deviation  and  loss  predictions  of  stator  blade  rows  are 
more  accurate  compared  to  the  prediction  of  losses  in  the  rotors.  In  all  cases  the  deviations  are  predicted  better  than  the  losses.  On 
the  whole  overall  mid-span  trends  at  almost  all  stations  are  predicted  well.  However,  there  are  discrepancies  between  measured 
and  calculated  values.  A  careful  examination  of  the  development  of  flow  through  the  blade  rows  shows  that  any  errors  in  total  loss 
and  deviation  distribution  in  the  spanwise  direction  propagate  downstream  and  can  be  detected  in  the  spanwise  property 
variations  at  the  downstream  stations. 
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5.  CONCLUSIONS 

The  goal  of  this  work  was  to  analyse  the  transonic  compressor  tests  available  in  the  open  literature  and  to  propose  possible 
improvements  in  total  loss  and  turning  correlations.  As  an  outcome  of  the  work  the  following  conclusions  are  reached: 

•  The  level  of  accuracy  of  deviation  prediction  should  be  greater  than  that  of  losses, 

•  The  off-design  deviation  correlation  which  gives  most  consistent  and  successful  results  was  found  to  be  that  of  Creveling. 

•  The  minimum  loss  incidence  predicted  by  the  NAS A-2D  correlation  needs  a  modification  separately  for  DCA  and  MCA 
blade  profiles.  This  modification  depends  on  the  blade  element  geometry  and  accounts  for  the  transonic  and  3-D  effects 
which  can  not  be  separated  in  the  present  investigation. 

•  The  design  deviation  angle  correlation  of  Carter  needs  a  correction  for  taking  into  account  the  transonic  and  3-D  effects 
which  are  not  separable  in  the  analysed  data. 

•  The  off-design  transonic  loss  correlation  developed  in  this  work  compares  quite  favourably  with  the  other  available  off- 
design  correlations  in  terms  of  overall  performance  prediction. 

•  Koch  and  Smith  [17]  minimum  loss  correlation  is  the  most  accurate  design  loss  prediction  method  available  in  the  open 
literature. 

The  comparison  of  the  experimental  results  of  a  2-stage  compressor  with  the  computed  results  using  the  new  correlation 
set,  revealed  that  the  new  set  can  be  regarded  as  satisfactory  in  giving  the  overall  trends  in  deviation,  losses  and  fluid  properties 
except  at  the  end  walls.  A  consistent  end-wall  boundary  layer  and  secondary  loss  calculation  method  is  necessary  for  more 
accurate  predictions.  The  loss  and  deviation  predictions  across  the  stators  are  more  accurate  compared  to  the  predictions  of 
losses  across  the  rotors.  For  more  accurate  predictions  of  the  flow  in  multi-stage  machines  spanwise  loss  mixing  procedures 
must  be  used.  Especially  for  off-design  predictions  the  minimum  loss  incidence  of  the  blade  row  together  with  the  flow 
direction  of  the  flow  ahead  of  the  blade  row  should  be  correct.  Small  errors  in  the  incidence  obviously  increases  in  a  cumulative 
way  through  the  compressor  in  the  case  of  multi-stage  machines. 
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APPENDIX  I 

NOTATION  AND  CASCADE  TERMINOLOGY 

A,  Annulus  area  of  the  streamtube 

a 

Ap  Annulus  area  contraction  ratio 

c  Chord  length 

Cq  Drag  coefficient 

C^  Skin  friction  coefficient 

Cp  Local  static  pressure  coefficient 

D  Diffusion  ratio 

h  Blade  height 

"k 

H  Form  factor  (6  /9) 

i  Angle  of  incidence  (°) 

La  Laval  number 

m  Correction  factor  for  deviation 

M  Mach  number 

r  Radius 

p  Pressure 

R  Gas  constant  for  air 

Re  Reynolds  number 

R^  Suction  surface  radius  of  curvature 

s  Pitch  length 

S  Entropy 

t  Blade  thickness 

V  Velocity  relative  to  blade  row 

V  Maximum  suction  surface  velocity 
max 

Meridional  velocity 
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Greek  Letters 

ot  Angle  of  attack  (°) 

B  Fluid  angle  (°) 

T  Ratio  of  specific  heats 

©  Deviation  angle  (°) 

★ 

6  Displacement  thickness  ,  Design  deviation 

r  Circulation  parameter 

e  Deflection  angle  (°) 

0  Momentum  thickness 

Prandtl-Meyer  expansion  angle  (°) 

K  Blade  angle  (°) 

V  Prandtl-Meyer  angle  (°) 

C  Stagger  angle  (°) 

p  Density 

a  Solidity 

$  Camber  angle  (°) 

u  Loss  coefficient 

0)  Rotational  speed  (rad/sec) 

Subscripts 

1  Inlet  condition 

2  Outlet  condition 

1s  Upstream  of  shock  wave 

2s  Downstream  of  shock  wave 

a  Sonic  condition 

c  Based  on  chord  length 

ch  Choke  condition 

cl  Clearance 

com  Compressible 

COR  Corrected  value 

CR  Critical  condition 
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eq  Equivalent 

inc  Incompressible 

LE  Leading  edge 

m  Average  (based  on  mean  velocity  vector  or  at  mean  radius) 

md  Modified 

min  Minimum 

p  Profile 

ref  Reference 

s  Shock 

St  Stall 

T  Total 

TE  Trailing  edge 

thr  Throat 

2  Axial  direction 

0  Tangential  direction 

Superscripts 

*  Design  condition 

"  Corrected  values 

Arithmetic  average 
Pitch  average 


A 
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BLADE  AND  CASCADE  TERMINOLOGY 


Fig. 1.1  Blade  Section  Profile  Terminology 

a  Location  of  maximum  camber 

d  Location  of  maximum  thickness 

c  Chord  Length 

Angle  between  blade-element  suction-surface  leading  edge  tangent 

line  on  conical  surface  and  meridional  plane  referenced  to 
leading  edge 

<25  Angle  between  blade-element  suction-surface  trailing  edge  tangent 

line  on  conical  surface  and  meridional  plane,  referred  to 

trailing  edge 

Angle  between  blade-element  mean  camber  line  on  the  conical 
surface  and  meridional  plane,  referred  to  leading  edge 

•"am  Angle  between  blade-element  mean  camber  line  on  the  conical 
surface  and  meridional  plane,  referred  to  trailing  edge 

$  Total  Camber  Angle 

tmax  Maximum  blade  thickness 

LER  Leading  edge  radius 

TER  Trailing  edge  radius 
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AXIAL 


Fig.  1.2  Cascade  Terminology 


s  Pitch  length 

Air  inlet  velocity 

'^2  Air  outlet  velocity 

Air  inlet  angle 
^2  Air  outlet  angle 

C  Stagger  (setting)  angle 

a  Solidity  (chord  length/pitch  length) 

E  Deflection,  difference  between  air  inlet  angle  and  air  outlet 

angle 
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i^,i  Mean  incidence  angle,  angle  between  inlet-air  direction  and  line 
tangent  to  blade  mean  camber  line  at  leading  edge 

6  Deviation  angle,  angle  between  exit-air  direction  and  tangent  to 

blade  mean  camber  line  at  trailing  edge 

The  dependent  parameters,  listed  above,  may  also  be  expressed  in 
arithmetical  form; 

®  "  '^1m  ■  '^2m 
c 

=  s 

e  =  -  ^2 

6  =  B2  -  <2m 

^m  =  ^1  -  ^1m 

Sometimes  it  is  convenient  to  define  incidence  angle  with  respect 

to  suction  surface.  As  shown  in  the  magnified  form  of  the  leading  edge 

(Figure  1.3).  Suction-surface  incidence  angle  is  denoted  as  the  angle 

between  inlet-air  direction  and  the  line  tangent  to  the  blade-suction- 

surface.  The  relation  between  i^^  and  i„  is  given  as 

ss  in  ^ 

^m  "  ^ss  *  ^"^IS  ■  "^Im^ 
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Fig.  1.3  Definition  of  Suction-Surface  Incidence  Angle 


47 


APPENDIX  II 

CORRELATIONS  AND  CALCULATION  PROCEDURE 


The  table  of  correlations  given  in  the  fol lowing, pages  are  compl ied 
from  the  relevant  references  on  design  and  off -design  loss  and  deviation 
correlations  present  in  the  open  literature.  The  first  column  of  the 
table  gives  the  references  from  which  the  correlations  are  taken.  The 
second  column  containes  the  equation  numbers.  Since  some  of  the  equations 
are  repeated  in  the  different  correlations,  such  equations  are  referred 
to  by  the  equation  numbers.  The  third  column  gives  the  correlation 
formulea  together  with  the  calculation  procedure.  The  sequence  of 
computation  is  indicated  by  giving  each  step  a  number.  In  most  cases, 
for  obtaining  the  required  result  from  a  correlation  several  empirical 
coefficients  are  needed.  These  coefficients  in  most  of  the  cases  are 
given  in  a  graphical  form.  The  figure  numbers  in  column  3  of  the  table 
are  the  figure  numbers  of  the  original  publications  (references  in 
column  1).  The  last  column  under  the  heading  "remarks"  include  informa¬ 
tion  with  regard  to  the  applicability  of  the  correlations  to  different 
flow  conditions  and  blade  geometry. 
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TABLE  OF  CORRELATIONS 


Ref. 


Eq*n 

No 


CORRELATION  FORMULAE  &  PROCEDURE 


REMARKS 


191 


(1) 


(2) 


PROFILE  LOSSES 


V2  |AV  I 

I)  Calculate,  D  «  I  -  +  -^gvY’ 


2)  »  f(D*)  (Fig. 148,  p.204)  and  H  =  1.1 


cos3i 


2H 

3H-I 


3)  w  =  2(^)  (cosB2^  '  ^  r.  .9^  oiTF^ 

^C-'  cosBaJ 


4)  Correct,  u)  =  f  (Re^»  Blade)  (Fig, 152,  p,206) 

u)  a  k  u)  ,  k  =  f  (M2,  H.  )  (Fig, 156) 
c  c  me  * 


Blading  NACA65  Series  Low 

speed  application. 

Correction  for  Re  &  M, 

Reasonable  results  up  to 

Re  <  2.5xl0^ 
c 

Design  only. 


1131 


(3) 


1)  Calculate, 


D 


eq 


1.12+a(a-a*)’^'^^+0.61 


cos^3i 


(tan3i-tan32) 


-I  cos32 
cos3i 


2)  (p  =  ^(^eq^  (Fig- 15)  and  H  =  1.08 

3)  Substitute  into  eq*n(2)  for  O) 


Incompressible , 
a  =  0.0117  (for  65-series) 
a  a  0.007  (for  C.4  series) 

Design  &  off-design. 


114] 


1)  (-)  ■= 


C  * 


1  -1  /  max  \ 


2)  The  values  of  k  &  C^*  from  Table  2,  H  =  1.08 

'  s  f 

3)  Substitute  into  eq*n(2)for  w 


Incompressible . 


V  occurs  near  the 
max 

leading  edge. 


116] 


1)  Calculate, 


(5) 


D 


COsB2rV  ,  I  -o  C0S^B2 

1.12+a(a-a*)  +0.61 


eq  cos3i 


riW 


r2 


Z2 
TzV 

.(tan3i-  — ^  tan32-  y —  (1- 

Zl  z 


o 


For  design  (a-a*)  =  0  and  D  =  D  * 

*  eq  eq 

2)  Correct,  (■^)  =  f  radial  pos.)  (Fig. 2) 


(6) 


3)  0)  =  2  (~). 


a  „  COs3l  2 
,  Q  ( _ \ 

C''*  COs32  *^COs32'' 


Off-design; 

(-)-(-)*“  (0.827Mi-2.692Mi^+2.675Mi^)(D  -D  *)%  D  >D  * 
c'^  ^c^  ^  eq  eq  ’  eq  eq 

(7)-(7)*=  (2.80Mi-8.71Mi^9.36Mi^)(D  -D  *)"  ,  D  <D  * 
c  c  eq  eq  eq  eq 


NACA65  and  C.4  series. 

Design  &  off-design.  For 
transoniccase  off-design 
depends  on  M.Spanwise 

0 

correction  for  (~)  Leakage  & 
secondary  losses  partly 
taken  into  account, 

a  a  0.017  (65  series) 

a  =  0,007  (C.4  series) 
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Ref. 


Eq’n 

No 


CORRELATION  FORMULAE  &  PROCEDURE 


REMARKS 


l29l 


(7) 


V2 

1)  Calculate.  ^  I- ^  (r.^-FOV.a 


2)  (|)*  =  f  (D+,  %  SPAN) 

(-)*  =  f  (D+,  %  SPAN) 
c 


for  rotors  (fig. 6) 
for  stators  (fig. 7) 


(8) 


3)  0)* 

"  in 


(7)* 


cos3 


Tip  clearance,  secondary, 
end-wall  B.L  losses  partly 
included. 


Desing  only. 


(9) 

(10) 


2 

1)  Calculate,  D*  =  1-  + 


.0* 


V.  AVg 
V  "  2aVi 


2)  ip  =  0.003  +  0.02375  D*  -  0.05  D*^  +  0.125  D*^ 

3)  Use  eq'n(6)for  tij 


ll8 


(11) 


Off-design  : 
1)  Determine 


(12) 


(13) 


(14) 


V  2 

^_maXx 
^  Vi  ^ 


where, 
V 


(1+  V  ^  ^ 


AV„ 


=  1+  (0.4+  |)  (^)+0.03+0.7  I 


2)  If  ,  correct  w  as, 

“cor  '  “  [ 


(15) 


3)  Correct  flow  angles  (if 

^Ist  =  ^ist  =  ^ich  =  ^Ich  "  1-5  ;  8.*"  =  Sa*  +  AS 


where,  A3  =  lO.Mi-7. 

4)  Then,  off-design  losses, 

w=  w*". (0.8333  +  0.1667  i  +  1.0) 

where , 

_  3i-3i*  _  3i  -3i* 

s  = -  for  3i  <  3i*  ;  s  = -  for  3i>3i* 


®lch-e>* 


NACA  65  series,  DCA,  PCA 
For  Transonic  Flows. 
Design  &  Off-design. 

MfCR  can  be  approximated 

^ICR  " 


For  off-design  operation, 
a  parabolic  variation  of 
loss  with  incidence  is 
assumed. 

Choke  and  stall  conditions 
are  reached  when,  oj  =  2.w 


(301 


(16) 

(17) 

(18) 


V  -V. 

1)  Calculate,  D*  =  - 

'1 

where, 

V  ,  AV„ 

max 


^.1.03+(0.4+-£)i-^+0.7f 


Incompressible . 


Blade  thickness  arbitrarily 
specified. 


2)  (!•)  -  (6.6  ^  +  0.34)(0. 0088+0. 0107D+-0.052D*“+0. 1160*^) 

3)  Use  eq*n(6)for 

4)  For  off-design,  follow  the  eq*ns.  (11),  (12),  (13),  (14) 
but  use  a  factor  of  1.8  instead  of  2.0  in  eq*n.(13) 


Design . 
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Ref. 


Eq*n 

No 


CORRELATION  FORMULAE  &  PROCEDURE 


REMARKS 


35 1 


(19) 


(20) 


For  Ml  <  1  , 

V  Vz 

1)  D  «  (-^)  -  *rr 

com  Vi  'com  Vi 

where, 


0  — ’J' 

2)  for  (— )  and  w  use  eq'ns.(18)&(6)respectively. 


For  Ml  >  1  , 


(21) 


(171 


1)  VJ' 


M,2 


(22) 


3)  Substitute  the  corrected  values  into  eq*ns.(17)&(20) 

V  „  V. 

4)  D 


^  maxv 

‘'com  ~  ^  Vi  'com 


5)  for  C-^)*  and  JI*  ,  use  eq* ns. (18)&( 6) respectively. 


l)  Calculate, 


(23) 


(24) 


(25) 


Vi  r 
D*  =  tj- 
eq  Vj  L 


C0s3i  2-11/0  ♦- 

(sln6i-Kiar*)2+( - )  ^^^(l+Ka  f +K,r*) 


where. 


r* 


02 


Ti+rj 

(—2—)  0  Vi 


A*  p* 
P  ^ 


M  ^  *  tan3i 

.  P*-  1- 

l-M  ^ 

Z1 


K2a7 


A  -A 


a;  -  (1-  ^>(1-  .  B = ^ 

cos  3  ai 


2)  D 


,♦  Fig^2a^  (1^  ^ 


eq 

♦ 

I 

eq 

4)  D*  ,  7-^  Fig.4a-b^  ^  (correction  for  streamtube 

^  height  variations) 

5)  If  Re^  >10®  ,  a  final  correction  for  Re^ 

~  Obtain  k  from  Appendix  2. 
s 

0 

-  Enter  Fig. 5  with  &  Re^  and  correct  (—) 

6)  Use  eq*n(2)for  calculating 


3)  D*  ,  Ml  >.  (~),  H  (correction  for  Mi) 


Design  only. 

Covers  up.  to  M^  *  1.4 

At  supersonic  inlet  flow, 

D  is  calculated  assuming 
an  equivalent  subsonic  flow 
ahead  of  the  shock  that 
leads  to  the  same  Mach 
number  at  the  position  just 
downstream  of  the  shock. 


Desing  only. 

Compressibility,  M,  Re, 
streamtube  contraction  are 
taken  into  account. 


Ki  »  0.2445 
Kz  =  0.4458 
Ka  =  0.7688 
U  *  0.6024 


0) 
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Ref. 

Eq'n 

No 

CORREUTION  FORMULAE  &  PROCEDURE 

REMARKS 

[20 1 

(26) 

(27) 

i-i* 

For  <  0 

acos^Bi 

^  ■  81+82 

cos^(  2  )  ^ 

A8*  -  ° 

acos^3i 

“t  ■  87+^ 

cos^(  2  ) 

0.017+0.91 

0.017+0.423  (^^)^‘^ 
.. 

5] 

Off-design  only. 

[19 1 

(28) 

(29) 

(30) 

1)  Calculate  w,j,  from  any  of  design  loss  correlations 

(the  method  described  in  reference  ll7l  is  recommended) 

2)  For  Ml  <  0.6 

Wj,  =  +  0.0005  (i-i*)^ 

3)  For  0.6  f  Ml  <  0.95 

ajj,  =  0^  +  (-0.0055+0.01  Mj)  (i-i*)* 

4)  For  Ml  >  0.95 

=  u*  +  (-0.0594+0.0667  Mj)  (i-i*)* 

Off-design. 

(31) 

(32) 

SHOCK  LOSSES 

^icR  =y s  y 

2)  (Cp,)^,„  =  1  -  (  “ 

where , 

1  -  (Cpi)  . 

axv 

j  '^inc 

Y-l/Y| ' 

[18] 

(33) 

V  AV^ 

(  r;)inc  =  (c  +  0-^)  °  Vi  +  0-7  c 

[17] 

(34) 

or, 

V  V  ^  AV. 

max  ~  max  _  i+o.7688  -  +  0.6024  0 
‘  thr  ^ 

(35) 

M  - 

yo.5[l+.-(.-l)La;,,]' 

[18] 

(36) 

Supercritical  subsonic  inlet  conditions 

Taking  into  account  the 
compressibility  effect  on 

0)  as  well  as  the  losses 

due  to  local  shocks. 

4)a)p  =  aj  +u)  =:(i)  .  (”2.0  (Mi-Mtp„)+i1 

'  r  p  s  p,inc  L  '  ICR^  J 
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Ref. 


Eq^n 

No 


CORRELATION  FORMULAE  &  PROCEDURE 


REMARKS 


[34] 


[11] 


(37) 


(38) 


(39) 


4)  0)™  =  0)  +0)  =  w  .  7i.4 

T  p  s  p,inc  ) 


Mi-(Micr-0.4)^ 


+1 


For  DCA  profiles  ®  .  is 

P,inc  ^ 

given  in  eq*n.8 


1)  M 


Is 


7  (M  +  1.0) 
2  Is 


2)  (1--^)  =  La 
^01 


1-^ 

2 

La¬ 

Y+1 

is 

1 

^  Y+1 

2 

La, 

Is 

1/Y-l 


The  shock  losses  are 
calculated  from  the 
pressure  distribution. 


(40) 

(41) 

(42) 


where  La  is  calculated  using  (31),  but  using 

the  values  at  shock  location  instead  of  minimum  ones, 


Lais  =  2  (L^ls+^) 


3)  Pitch  averaged  total  pressure  at  the  exit  of  the  cascade, 


^02  *^s  ^  /  ^02s'\ 

Pqi  =  -M,  j 


[16] 


Supersonic  inlet  conditions; 


(43) 


(44) 


1  ^  D  le 

1 )  Ru  =  — + 


^raax  ^LE  fc  ^LE-i  .  rc  ^LEi  ^ 

1 - 2~  L  2  -  ~J  ^4^^  ^  2  -  ^ 


,3n(|)  . 


6  ^  .  *'max  *'LE 


2)  The  Prandtl-Meyer  expansion  angle  is, 
s.  sin3i 


QpM  =  tan 


-1 


5  3i+Ru 


(45) 


(46) 


(47) 


(48) 


(49) 


3)  The  Prandtl-Meyer  angle  is, 

Vi  =  ^  {tan“^  ^  (Mi"-!)^''^  }  -  tan“^(Mi "-1)^''^ 

4)  After  the  wave,  the  Mach  number  is  given  by, 

«ls  =  =  ^(^s) 

or. 


M^g  =  1.0  +  (0.0432) 


5  )  The  normal  shock  entry  Mach  number  is. 


«ls  =  I  («^+«ls) 


6  )  From  normal  shock  tables,  with  M.  ,  find  the  value 


p  -P  1-p  /p 

01s  02s  02s'  01s 


^Ols-^^ 


1-Pi  /P, 


01s 
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Ref. 

Eq*n 

No 

CORRELATION  FORMULA  &  PROCEDURE 

REMARKS 

[17] 

(50) 

1 )  Losses  due  to  leading  edge  bluntness  is  given  in  terms  of 
entropy  rise, 

-In  {l-tj^j,/(s.cosBi).[l.28(Mi-l)+0.96(Mi-l)^]} 

2 )  Shock  losses  due  to  normal  and  oblique  shocks  are  given  in 
figure  7. 

There  exist  3  types  of 

shock  losses, 

1-  Due  to  leading  edge 

bluntness 

2-  Due  to  normal  passage 

shock 

3-  Due  to  oblique  shock. 

(32) 

(51) 

DEVIATION  , 

1 

1)  Base  equation  is  given  as 

6*  s  ni4)  /  — 
c 

2)  m  =  ra  (profile  shape,  C  )  Fig,  A, 3  in  [11]  or  Ref.  [32] 

(291 

(52) 

1)  Modified  base  equation  is  given  as 

(e-i)  m  , 

X*  ^  ^  md  c 

md  y — 

1  -  m  , 

rad  c 

a  2  - 

2)  m  ^  =  0.92  (-)  +  0.002  02 

^  md  c 

or, 

e  fluid  turning  angle 

degrees, 

02  average  fluid  angle  at 

cascade  exit  measured  from 

axial  direction,  degrees. 

[33] 

(53) 

=  (0.219  +  0.0008916  S  +  0.00002708  5^). 

175-0. 03552C  +  0.00019175^ 
c' 

[9] 

(54) 

(55) 

(56) 

1)  Base  equation  is  given  as, 

6°  ,  =  6°  +  (t> 

ref  0  ^6 

2)  «o  =  (''6>sh 

3)  For  off-design 

6°  =  6°  ,  +  (i-i  ,)(^)  , 

ref  ^  ref  ^di'ref 

(6°)io=  f(B.o)(Fig.l61,SP-36) 
”0=1=  (Fig.l63,SP-36) 

b  =  f(Bi)  (Fig. 164,  SP-36) 

(K5)^=  f(^)  (Fig. 172,  SP-36) 

(K^)sh=  f (Blade  shape) (SP-36) 
=  1,1  (for  C  series) 

(^)ref=f(Bi.a)(Fig. 177, SP-36) 

[11  ] 

(57) 

(58) 

(59) 

(60) 

(61) 

(62) 

1)  Base  equation  is  given  as 

f  K3(t>  1 

e  -  K1.K2  [  +  BJ 

where, 

K3=  (5-2  !■)  / 5^-1000  +  100  (5.5  -  2.6  ^) 

B  .  8  (-)^  -  17  -  +  16 
c  c 

Ki  -  1-0.28  (x^  -  0.40)  for  -  0.3  to  1.0 

K2  =  1-0.016  (10-c)  for  c  =  1.25  to  12.5% 

at  base  incidence 

.  -1 

1  =  Ki  «  -  Sin  — 

0  ^m,0  s 

-  Maximum  thickness 

Length  of  camberline 

-  Location  of  max, thickness 
c  Lenth  of  camberline 

0*  Blade  angle  at  LE 
*  (measured  from 

tangential  direction) 

Fp  *  minimum  flow  blade 

passage  width  (throat) 
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Ref. 

Eq'n 

No 

CORRELATION  FORMULAE  &  PROCEDURE 

REMARKS 

l)  Correct  design  air  inlet  angle  for  Mach  number, 

For  off-design. 

(63) 

(64) 

(Bi*)'  =  Bi*  +  10  Ml  -  7.0  for  Ml  >  0.7 

(Bi*)'  =  Bi*  for  Ml  <  0.7 

The  method  suggested  by 
Lieblein  1 13l  is  used  with 
corrections  applied  for 

2)  Corrections  for  streamline  slope  &  axial  velocity  variation 

Ml  ,  streamline  slope  and 

are  applied  as, 

axial  velocity. 

(65) 

2  tan  ["(Si*)  Vcosal 

•^«"(e?^or)  =  W 

(1+/^) 

mi 

a  =  streamline  slope. 

(18) 

3)  Assuming  the  air  turning  angle  varies  linearly  with  incidence 

in  minimum-loss  region, 

(66) 

{AS*  +  (i-i*)) 

where,  A6*  =  (Ji  -  6*  +  i* 

(67) 

rdX|B11  *  ^  0.4805  +  (0. 00407-1. 38xlO“^)Bi* 

L  di  J  COR 

+  (0.698-0.005  B*  +  1.51x10"^  B*^  )a 
*C0R  COR 

+  (-0.226+0.00152  B*  -  0.431x10"^  6*^)a^ 

*C0R  COR 

(68) 

4)  tan32  -  2^W  /W  +1^ 
m2  mi 

(69) 

5)  6  =  B2-  C^i-W 

[161 

(70) 

1).  6  =  6*+  [6.40-9.45  (Mi-0.6)]  (Deq-D*q) 

For  off-design 

For  DCA  Blading  6*  may  be 
obtained  from  either  i9l 
or  [32I 

(71) 

6-6* 

1)  -  f(x) 

6B 

2)  for  X  >  0 

For  off-design  . 

(72) 

f(x)  =  -0.809xl0"^+0.5588x-0.2928x^ 

[19] 

for  X  <  0 

(73) 

f(x)  =  0.1191xlO~\o.480x+0.3452x^ 

where , 

(74) 

i-i* 

X  «  * 

AB* 

(75) 

1)  32*  S1-A6*  {f(x)} 

2)  for  x  >  0  , 

(76) 

f(x)  =  1.0  +  0.86x  -  1.36  x^ 

For  off-design. 

[20] 

for  X  <  0  , 

(77) 

f(x)  =  1.0  +  0.86x  -  0.35  x^ 

3)  where  x  is  defined  by  equation  (74) 
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APPENDIX  III 

DATA  UNCERTAINTY 

Table  III.1  lists  the  precision  errors  given  in  the  reports  on 
which  the  work  explained  in  this  report  is  based  on.  It  was  not  possible 
to  apply  uncertainty  analysis  to  the  averaging  procedures  due  to  the 
lack  of  relevant  data.  However,  the  uncertainty  limits  given  in  the 
table  is  propagated  through  the  equation  2.2,  (Section  2)  which  is  used 
for  loss  coefficient  calculations.  The  precision  index  for  the  loss 
coefficient  (S^  )  is  calculated  using  Taylor  series  expansion  [43] 


The  partial  derivatives  in  the  above  equation  are  evaluated  from 
equation  2.2  and  precision  index  of  the  loss  coefficient  is  calculated 
for  197  sample  points  of  reference  [1].  If  was  found  out  that  the  maximum 
occurance  is  at  the  precision  error  of  +  1.78  %  on  the  total  loss.  In 
this  calculation  the  precision  error  of  r^  and  are  taken  as  +  0,33 
mm. 


For  analysing  the  effect  of  mislocating  the  streamline,  a  single 
operating  point  is  chosen  and  it  is  assumed  that  the  upstream  location 
of  the  streamline  is  correct.  The  downstream  radius  is  then  changed 
between  +  0.5  cm  along  the  span.  For  each  downstream  spanwise  streamline 
location,  thermodynamic  quantities  measured  at  the  corresponding 
location  are  taken.  The  precision  errors  on  the  measured  quantities  and 
on  r 2  are  then  used  in  Equation  111.1.  The  change  of  uncertainty  interval 
due  to  the  mislocation  of  the  downstream  point  of  the  streamline  can  be 
seen  in  Figure  111.1,  The  figure  shows  that  the  precision  index  S_ 
reduces  +  0,04  %  when  the  streamline  position  is  misplaced  by  0.5  cm 
towards  the  bub  of  the  blade.  If  the  streamline  is  misplaced  by  0,5  cm 
towards  the  tip  of  the  blade  the  precision  index  on  Sj  increases  by 
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APPENDIX  IV 

SUBROUTINE  FOR  THE  NEW  LOSS  AND  DEVIATION  SET 


The  following  subroutine  is  prepared  in  FORTRAN  IV  language  and 
used  in  the  throughflow  computer  code  for  performance  prediction. 
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SUBROUTINE  CASCADE! AAl ,  AA2,  ALFA, CAM, CODE,  CORD, ETl, ET2, 
lOM,  H1H2,  LER,  MIN,  OMEGA,  PHI,  PH2,  PI,  POl,  Rl,  R2,  ROl,  SOL,  TlX, 
2TRAT,  WMl ,  WM2,  WZ 1 ,  WZ2,  XBETl ,  XKAP 1 , DEV,  PLOD > 

IMPLICIT  REALMS < A-H, 0-Z) 

DIMENSION  H0HM(4),  T0TM(4),  DFACT(4), H0HH(4> 

REAL*8  MIN,  KAPl,  INCH 


C-- 

c 

n 

FOLLOWING  DATA  SHOULD  BE  SUPPLIED  TO  SUBROUTINE  CASCADE 

Lf 

c 

AAl 

ANNULUS  AREA  AT  THE  INLET 

c 

AA2 

ANNULUS  AREA  AT  THE  EXIT 

c 

ALPA 

STREAMLINE  SLOPE 

c 

CAM 

CAMBER  ANGLE 

c 

CODE 

1  SPECIFIES  ROTORS  ;  2  SPECIFIES  STATORS 

c 

CORD 

CHORD  LENGTH 

c 

ETl 

DIMENSIONLESS  RADIUS  AT  THE  INLET 

c 

ET2 

DIMENSIONLESS  RADIUS  AT  THE  EXIT 

c 

GM 

RATIO  OF  SPECIFIC  HEATS 

c 

H1H2 

STREAMTUBE  HEIGHT  RATIO 

c 

LER 

BLADE  LEADING-EDGE  RADIUS 

c 

MIN 

RELATIVE  INLET  MACH  NUMBER 

c 

OMEGA 

ANGULAR  SPEED 

c 

PHI 

FLOW  COEFFICIENT  AT  THE  INLET 

c 

PH2 

FLOW  COEFFICIENT  AT  THE  EXIT 

c 

PI 

INLET  STATIC  PRESSURE 

c 

POl 

INLET  TOTAL  PRESSURE 

c 

Rl 

RADIUS  AT  THE  INLET 

c 

R2 

RADIUS  AT  THE  EXIT 

c 

ROl 

INLET  DENSITY 

c 

SOL 

SOLIDITY 

c 

TlX 

INLET  TEMPERATURE 

c 

TRAT 

BLADE  THICKNESS  TO  CHORD  RATIO 

c 

WMl 

MERIDIONAL  VELOCITY  AT  THE  INLET 

c 

WM2 

MERIDIONAL  VELOCITY  AT  THE  EXIT 

c 

WZl 

AXIAL  VELOCITY  AT  THE  INLET 

c 

WZ2 

AXIAL  VELOCITY  AT  THE  EXIT 

c 

XBETl 

AIR  ANGLE  AT  THE  INLET 

c 

c 

c — 

XKAPl 

BLADE  ANGLE  AT  THE  INLET 

c 

c 

OUTPUTS  ARE  ; 

c 

DEV 

:  OFF-DESIGN  DEVIATION  ANGLE 

c 

c 

r — 

PLOD 

:  OFF-DESIGN  LOSS  COEFFICIENT 

c 

c 

KOUT  SPECIFIES  OFF-DESIGN  DEVIATION  ANGLE  CALCULATION 

METHOD 

c 

K0UT=>1 

MOFFAT  <  L I EBLE I N+C OR  REC  T I ONS ) 

c 

K0UT«2 

SWAN 

c 

K0UT«3 

C REVELING 

c 

c 

c — 

K0UT«4 

HOWELL 

c 

c 

KSER  SPECIFIES  THE  TYPE  OF  PROFILE 

c 

KSER=1 

DC A  PROFILES 

c 

KSER =2 

MCA  PROFILES 

C 
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PRELIMINARIES 

C  a- »»»»##**»»»»»»»•»••»••»••»••»••»••»••»•»••»••»••»••»••»••»••»••»••»••»••»••»••»••»••»••»••»•*»•»•■«•»•  a  •»•■» -a -» -s  •»•»•«••«••«••»••»•♦•♦•»  •» 

C 

PI=3. 141592 
DEG=PI/180. 

INCH»100.  /2.  54 

K0UT=2 

KSER=2 

SH=0.  7 

EQ=0.  007 

BET1=DABS<XBET1 > 

KAP1=DABS(XKAP1 > 

STG=DABS(KAPl-CAM/2.  > 

ACIN=XBET1-XKAP1 
XKS“0. OOOlaiNCH 
C0RD=CaRD4HNCH 
RTIP1=R1/ET1 
RTIP2*R2/ET2 
C 

Caaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa*  aaa 

Caaaaaaaaaa'^aaaaaaaaaaa  DESIGN  INCIDENCE - AINaaaaaaaaaaaaaaaaaaaaaa* 

Caaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa-* 

C 

C -  SP-36  DESIGN  INCIDENCE  (MODIFIED) 

B1M=»KAP1 

SLP»<-.  063+.  02274aS0L)  +  (-.  0035+.  0029aS0L )aBlM+ ( -.  379E-4-1.  llE-5* 
lS0L)aBlMa*2 

ZCI*(. 0325-.  0674aSaL)  +  (-.  2364E-2+.  0913aS0L)aBlM 
1+(1. 64E-5-2.  38E-4aS0L)aBlM»a2 
C0R=-6.  19E-3+13. 64aTRAT-93.  7aTRATaa2+274.  3aTRATa*3 
AIN=CORaSH'aZCI+SLPaCAM 
GO  TO  (100»200)  KSER 
100  AIN=AIN+0.  7238aMIN+7.  5481 
GO  TO  300 

200  AIN=AIN+1.  3016aMIN+5.  7380 
300  B1M=KAP1+AIN 
C 

Caaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa-»--»-»--»-aaaaaaaaaaaaaaxaaaaaaaaaaaaaaaa-*a 

Caaaaaaaaaaaaaaaaaaaaa  DESIGN  DEVIATION - DEVM  aaaaaaaaaaaaaaaaa**#*-*-* 

Caaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa  a  aaaaaaaaaaaa* 
C 

Q -  CARTERS  DEVIATION  RULE  (MODIFIED) 

SC=0.  216+9.  72E-4aSTG+2.  38E-5*STG«a2 
DEVM=SCaCAM/DSQRT(SOL) 

DEVM«-.  1988183*DEVMaa2+3.  0186aDEVM-l.  099379 
C 

Caaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa  aaaaaa  aaaaaaaaaaaaaa 
Caaaaaaaaaaaaaaaaaaaaaaaaaa  PRELIMINARIES  aaaa  >aaaaaaaaaaaaaaaa#*»aaa-» 
Caaaaaaaaaaaa*aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa»>>aa»»aaa»aaaa*»ii'a 
C 

C  - COMPUTE  THE  FOLLOWING  DESIGN  PARAMETERS - 

C  Blai  B2a.  DELBa.  DEQa 

DELBM=C AM-DE VM+ A I N 
B2M=B1M-DELBM 
BlR=BlMaDEG 
B2R=B2MaDEG 
Bl=BETlaDEG 
CSB1=DC0S(B1R) 

TNB1*DTAN(B1R) 


a  a 
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CSB2=DC0S(B2R) 

TNB2=DTAN(B2R) 

C  - MINIMUM  LOSS  DIFFUSION  FACTORS - 

AK1=ET2»PH2/(ET1»PH1 ) 

C  - CORIOLLIS  TERM - 

AK2=ET1 /PHI* (1.0- {ET2/ET1 ) **2 ) 

IF(CODE.  EQ.  2)  AK2=0.  0 
AK3=EQ*( DABS (AC IN-AIN) )*1.  43 
V1V2=PH1*CSB2/(PH2*CSB1 ) 

C 

C*************************************-****#**!! 

C************  DESIGN  TOTAL  PRESSURE  LOSS  COEFFIC IENT**^J **>*********» 

C*****************************»»»»»»»»»»#»»»»»-»»>J|.Am.^t^f»»».»#.!HJ.»#»**4J.* 


- KOCH  y.  SMITHS  CORRELATION  TO  COMPUTE  DESIGN  LOSS - 

IF(CODE.  EQ.  2)  0MEGA=0.  0 
VTET2=0MEGA*R2-WM2*DTAN(B2R ) 

VTET1»0MEGA*R1-WM1*DTAN(B1R ) 

V1=DSQRT(VTET1**2+WM1**2) 

IF(CODE.  EQ.  1)  V1=DSQRT( (WM1*DTAN(B1R) )**2+WMl»*2) 

C  - CALCULATE  BLADE  CIRCULATION  PARAMETER - 

RAVE*(Rl+R2)/2 

C IRC= (R 1 *VTET1-R2*VTET2 ) / ( SOL*R AVE*V1 ) 

BETAV»(BlM+B2M)/2.  0 
BETAVR=BETAV*DEG 
XMZ1=WZ1/DSQRT(GM*287.  05*T1X) 

C  - CALCULATE  ANNULUS  AREA  CONTRACTION  RATIO - 

AACR=(1.  0-0.  4438*S0L*TRAT/DC0S(BETAVR) )*( 1.  0- ( ( AA1-AA2 ) / 

1 (3. 0*AA1 ) ) ) 

R0PR01  =  1.  -(XMZ1**2/(1.  -XMZ1**2) )*(1.  -AACR-.  2445»(DTAN(B1R) / 
IDCOS(BIR) )»SOL*CIRC) 

VMAXVP=1.  0+0.  7688*TRAT+0.  6024*CIRC 

VPVl»((DSIN(BlR)-0.  2445*S0L*CIRC)**2+(DC0S(B1R)/(AACR»R0PR01 ) )**2 
+  )**0.  5 

C  - CALCULATE  DESIGN  DIFFUSION  FACTOR - 

DEQM*VPV1*VMAXVP*V1V2 
IF(DEQM.  GT.  2.  )  DEQM=2.  O 

C  - CALCULATE  DESIGN  MOMENTUM  THICKNESS - 

DMTU=0.  0216*DEQM**2-0.  0557*DEQM+0.  0406+0.  0025 

C  - CALCULATE  FORM  FACTOR - 

HTEU*2.  047628*DEQM**2-4.  902983*DEQM+4.  2060613 
IF(DEQM.  GE.  2.  )  HTEU=2.  6 

C - CORRECT  DMT  V.  HTE  FOR  INLET  MACH  NUMBER - 

HOHMd  )=0.  297619*MIN**2+0.  1 09226*MIN+0.  99 
H0HM(2)=0.  188988*MIN**2+0.  069702*MIN+0.  994 
H0HM(3)=0.  140557*MIN**2+0.  039275*MIN+0.  9992 
H0HM(4)=0.  088474*MIN**2+0.  035108*MIN+0  9984 
C 

T0TM(l)=-0.  045103*MIN**2-0.  018419*MIN+1 .  0035 
T0TM(2)«-0.  08846*MIN**2-0.  01967*MIN+1.  004 
T0TM(3)=-0.  1 10173*MIN**2-0.  033056*MIN+1 .  004 
T0TM(4)=-0.  120996*MIN**2-0.  065324*MIN+1  006 
C 

DFACT(1)  =  1.  0 
DFACT(2)=1.  3 
DFACT(3)=1. 5 
DFACT(4)=1. 7 
C 

CALL  INTERP ( DFACT#  HOHM<  DEQM, HOHMX#  4) 

CALL  INTERP (DFACT. TOTM<  DEQM, TOTMX. 4) 
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C 

IF(DEQM.  GT.  1.  7)  H0HMX=1.  0 

HTE»HTEU»HOHMX 

DMT=DMTU*TaTMX 

C - CORRECT  DMT  J»y.  HTE  FOR  STREAMTUBE  C0NTRACTI0N(H1H2) - 

DMT«DMT»(0.  588235»HlH2+0.  411765) 

H0HH(l)*-0.  02143»H1H2+1.  02143 
H0HH(2)=-0.  00714»H1H2+1.  00714 
H0HH(3)=0.  04524»HlH2+0.  75476 
H0HH(4)=0.  157l4*HlH2+0. 84286 
CALL  INTERPCDFACT.HOHH, DEQM, HOHHX, 4) 

HTE»HTE*HOHHX 

C  - CORRECT  DMT  FOR  REYNOLDS  NUMBER - 

XMU1*0.  00001*(0.  0045833»TlX+0.  468751 ) 

XNU1»XMU1/R01 
REX=(V1»C0RD/INCH)/XNU1 
RROUGH=XKS/CORD 
IF(RROUGH.  NE.  0.  )  GO  TO  222 
548  IF(REX.  LE.  200000.  )  T0TR=589.  6»REX»* (-0.  5 ) 

IF(REX.  GT.  200000.  )  T0TR=10.  0»REX»»(-0.  166) 

GO  TO  333 
222  CONTINUE 

IF(RROUGH.  GE.  0.  005)  GO  TO  111 
IF(RROUGH.  GE. 0. 001 )  GO  TO  544 
IF(RROUGH.  GE.  0.  0005)  GO  TO  545 
IF(RROUGH.  GE.  0.  0001 )  GO  TO  546 
IF ( RROUGH.  GE. 0.  00005)  GO  TO  547 
547  THTH»( (RROUGH-0.  001 )/0.  004)*0.  115+0.915 
IF(REX.  LT.  1800000.  )  GO  TO  548 
TOTR=THTH 
GO  TO  333 

546  THTH=( (RROUGH-0.  001 )/0.  004 )*0.  36+1.  03 
IF(REX.  LT.  900000.  )  GO  TO  548 
TOTR=THTH 
GO  TO  333 

545  THTH»( (RROUGH-0. 001 )/0. 004)*0.  58+1.  39 
IF(REX.  LT.  108000.  )  GO  TO  548 
TOTR=THTH 
GO  TO  333 

544  THTH=((RR0UGH-0..001)/0.  004) *2.  33+1.  97 

IF ( REX.  LT.  90000.  )  GO  TO  548 
TOTR=THTH 
GO  TO  333 

111  IF(REX.  LT.  18000.  )  T0TR=10.  0»REX*»  (-.  166) 

IF(REX.  GE.  18000.  )  T0TR=4.  3 
333  CONTINUE 

C  - CORRECT  HTE  FOR  REYNOLDS  NUMBER - 

H0HR=2.  29087»REX»»(-0.  06) 

DMT=DMT»TOTR 
HTE*»HTE*HOHR 
C  - CALCULATE  PLM - 

AA=( (2.  »HTE)/(3.  »HTE-1.  ))/((!.  -DMT*SOL*HTE/DCOS( B2R ) )**3) 
PLM=2.  »DMT*(S0L/DC0S(B2R)  )»(  (DC0S(B1R)/DC0S(B2R)  )«»2.)»AA 
C  - CALCULATE  SHOCK  LOSS - 

PLS=(P01/(P01-P1 ) )»( (2»S0L»LER/C0RD»DC0S(BIM) )#(1. 28(MIN- 
11.  )+0.  96(MIN-1.  )»»2)  ) 

PLM=PLM+PLS 

C  CORRECT  Bl*  FOR  MACH  NO 

IF(MIN.  LE.  0.  7)  GO  TO  10 
DB=10.  »MIN-7. 
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Biri=BIf1+DB 
10  CONTINUE 
C 

OFF-DESIGN  DEVIATION 
C 

GO  TO  ( 483,  486,  487,  488 ) ,  KOUT 

485  CONTINUE 
C 

C  - LIEBLEINS  OFF-DESIGN  DEVIATION  CORRELATION - 

C  - CORRECT  FOR  AXIAL  VELOCITY  VARIATION  AND  SLOPEdlOFF AT) - 

B 1 E=B IM/DCOS ( ALP A»DEG ) 

VRAT=WM2/WM1 

B1E=DATAN(2.  »DTAN(B1E»DEG> / ( 1.  +VRAT) ) /DEG 

DBDI*0.  4805+(0.  00407-0.  000138»B1E  )»BlE+(0.  678+(-0.  005 >• 

10.  000151»BlE)»BlE)»S0L+(-0.  226+(0.  00152-0.  431E-4*B1E)  »B1E)«SQL^2 
Z I =AC IN-AIN 

B2E=BETl-(DELBri+DBDI»ZI ) 

BET2=DATAN(0.  5*(1. +1.  /VRAT)*DTAN(B2E»DE0*DC0S(ALPA*DEG) ) )/DEG 
GO  TO  487 

486  CONTINUE 
C 

C  - SWANS  OFF-DESIGN  DEVIATION  CORRELATION - 

BET2=DEVM+ ( KAP 1 -C AM ) 

KNT=1 

250  CONTINUE 

B2=BET2*DE0 

VRM=1.  12+AK3+0.  61»DC0S( B1 ) *»2» ( DTAN (B 1 )-AKl*DTAN( B2)-AK2 > /SDL 
V1V2=PH1»DC0S(B2)/(PH2»DC0S(B1> ) 

DEQ»VRM*V1V2 

C 

C - IF  FLOW  SEPARATES,  USE  A  LIMITING  VALUE  OF  DEG - 

IF(DEQ.  GT.  2.  2)  DEQ=2.  2 

DIFF=DEQ-DEQM 

AM=MIN 

IF(AM.  LT.  0.  65)  AM=0.  65 
DF=(6.  4-7.  45»(AM-0.  6)  )»DIFF 
DEV*DF+DEVM 
BET2=DEV+(KAP1-CAM) 

IF(KNT.  EQ.  2)  GO  TO  260 
KNT=KNT+1 
B2=BET2»DEG 
GO  TO  250 
260  CONTINUE 
GO  TO  487 

487  CONTINUE 
C 

C  - CREVELINGS  OFF-DESIGN  DEVIATION  CORRELATION - 

X= (AC  IN-AIN )/DELBM 

IF(X.  GE.  0.  0)  F=0.  807E-3+(.  5588-.  2728»X)»X 
IF(X.  LT.  0.  0)  F=0.  1171E-3+(.  48+0.  3452»X)»X 
DEV=DELBM*F+DEVM 
BET2=DEV+KAP 1-C AM 
GO  TO  487 

488  CONTINUE 
C 

C  - HOWELLS  OFF-DESIGN  DEVIATION  CORRELATION  — 

X=(AC IN-AIN) /DELBM 
Fl  =  l.  36 


n  n  n  n  n  n 
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IFCX.LT.  0.  0)  F1=0.  35 
F=l.  0+0.  86*X-F1*X**2 
BET2«-DELBM*F+BET1 
DEV«BET2-KAP 1+C AM 
489  CONTINUE 

C0B1»DC0S(B1 ) 

TABl^DTANCBl ) 

B2=BET2*DE0 

C0B2=DC0S(B2> 

TAB2=DTAN<B2> 

DEV=BET2- <  KAP 1 -C AM ) 

C 

c*****************  OFF-DESIGN  LOSS  COEFFICIENT  *•«*»*#•»**■»» 

C*****************'»nnnnnnnnnnnnnnf»********»##5i  >  »•»*•»##*««#•!»■♦ 


- THE  NEW  OFF-DESIGN  LOSS  CORRELATION - 

AX=AC IN-AIN 

GO  TO  <450. 500)  KSER 

450  IF(MIN.  LT.  0.  56)  GO  TO  451 
IF(AX.LT.  0.  )  AEX=0.  05336*MIN-0.  02937 
GO  TO  452 

451  AEX=0. 

452  IF(AX.GE.  0.  )  AEX=0.  00500#MIN-0.  00075 
GO  TO  550 

500  IF<MIN.  LT.  0.  62)  GO  TO  501 
IFIAX.LT.  0.  )  AEX=0.  02845»M IN-0.  01741 
GO  TO  502 

501  AEX=0. 

502  IF(AX.GE.  0.  )  AEX=0.  00363*MIN-0.  00065 
550  PL0D»AEX»AX**2+PLM 

RETURN 

END 

-  SUBROUTINE  INTERP  INTERPOLATES  THE  PARAMETERS  - 

SUBROUTINE  INTERP (PS I, PROP, PSIX. PROPX, NYY) 

IMPLICIT  REALMS ( A-H, O-Z) 

DIMENSION  PSKNYY).  PROPINYY) 

NYNY*NYY-1 

1  IFIPSIX.  LT.  PSHD)  GO  TO  3 
IF(PSIX.  GT.  PSKNYY)  )  GO  TO  4 
DO  5  I  =  1,NYNY 

IF(PSIX.  GE.  PSK  I  ).  AND.  PSIX.  LE.  PSI  ( I  +  l )  )  GO  TO  6 

5  CONTINUE 

3  PR0PX*PR0P<1 ) 

RETURN 

4  PROP X»PROP< NYY) 

RETURN 

6  PROPX«(PSIX-PSI  ( I+l )  )/<PSI  ( I  )-PSI  ( I  +  l )  )'»PROP<  I ) 
PROPX*PROPX+<PSIX-PSI ( I ) )/(PSI ( I+l )-PSI < I ) )»PROP< I+l ) 
RETURN 

END 
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